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SUMMARY 


Stresses  were  calculated  for  finite-size  orthotropic  laminates  loaded  by 
a  frictionless  steel  pin  in  a  circular  hole  of  the  same  diameter.  These  cal-  / 
culations  were  based  on  finite-element  analyses  for  six  laminates;, j  quasi- 
isotropic  [0°/±45°/90O] g,  0°'  90°,  [0°/90O]s,  [±450]s,  [0O/±45°]s.  Stress 

concentration  factors,  based  on  nominal  bearing  stress,  were  determined  for  wide 
ranges  of  the  ratios  of  width  to  diameter  and  edge  distance  to  diameter 

(^dl)  For  reference  purposes,  the  infinite-laminate  case  was  also  analyzed  for 
each  laminate.^ 


Orthotropy  had  a  significant  influence  on  both  the  magnitude  and  location 
of  the  maximum  tensile  stress  concentration  on  the  boundary  of  the  hole.  For 
example,  the  stress-concentration  factors  for  the  infinite-laminate  cases  ranged 
from  0.82  to  2.16,  compared  with  0.87  for  the  quasi- isotropic  laminate.  The 
laminates  with  0^  plies  developed  the  peak  tensile  stress  near  the  ends  of  the 
pin/hole  contact  arc.  But  the  90^  and  ±45^  laminates  had  peaks  where  ply 
fibers  were  tangent  to  the  hole. 

|The  finite  widths  and  edge  distances  strongly  influenced  the  tensile 
stress  concentration.^  For  the  practical  range  w/d  ^  3,  the  peak  tensile 
stresses  were  as  inuch'as  80  percent  larger  than  the  inf inite- laminate  reference 
value.  Even  for  rather  wide  laminates  with  w/d  =  10,  some  of  the  peak  tensile 
stresses  were  amplified  by  nearly  20  percent.  For  e/d  ^  3,  these  stresses 
were  amplified  by  as  much  as  50  percent.  In  contrast,  the  finite  width  and 
edge  distance  had  little  effect  on  shear-out  and  bearing  stress  concentrations. 
For  w/d  ^  3  and  e/d  ^  3,  the  peak  shear-out  stress  was  amplified  by  less 
than  25  percent  and  peak  bearing  stress  by  only  10  percent. 


INTRODUCTION 

The  stress  concentrations  at  fastener  holes  are  widely  recognized  as  having 
a  degrading  effect  on  structural  fatigue  life.  In  conventional  metallic  struc¬ 
tures,  fatigue  cracks  commonly  initiate  at  these  stress-concentration  sites. 

In  recognition  of  this,  the  effects  of  stress  concentrations  on  crack  initiation 
and  growth  have  been  incorporated  into  current  methods  for  predicting  fatigue 
life.  However,  in  composite  structures,  the  stress  concentration  effects  are 
much  more  complicated  and  are  not  well  understood.  For  example,  holes  can  pro¬ 
duce  either  composite  delamination,  matrix  crazing,  fiber  failures,  or  various 
combinations  of  these.  Furthermore,  for  a  pin-loaded  hole,  these  different 
types  of  damage  can  produce  three  different  modes  of  failure:  net- tens ion, 
shear-out,  or  bearing.  Before  analytical  models  can  be  developed  for  these 
various  failures,  the  stress  concentrations  at  loaded  holes  must  be  better 
understood.  To  date,  the  few  stress  analyses  that  have  dealt  with  this  problem 
have  been  limited  to  rather  specific  cases.  The  purpose  of  the  present  paper 
is  to  present  some  of  the  needed  background  on  stresses  at  loaded  holes  in 
compos i te  laminates. 


This  study  was  based  on  a  two-dimensional  finite-element  analysis  of 
stresses  in  an  orthotropic  laminate  that  is  loaded  by  a  frictionless  steel  pin 
in  a  circular  hole  of  the  same  diameter »  Six  graphite/epoxy  laminates  were 
analyzed:  quasi- isotropic  [o°/±45^/90^ 0^,  90^,  [o^/90^]s,  [+45^  ]g,  and 

[o^/±45'^]3.  Finite-element  models  were  developed  for  wide  ranges  of  width- to- 
diameter  ratio  v;/d  and  edge-distance-to-diame ter  ratio  e/d.  Firsts  each 
laminate  was  analyzed  for  different  values  of  vz/d  using  one  large  value  of 
e/d.  Next,  a  range  of  e/d  was  used  with  one  large  w/d  value. 

The  results  are  presented  in  two  ways:  first  as  stress  distributions  near 
the  hole  and  then  as  stress-concentration  factors  for  ranges  of  w/d  and  e/d. 
Results  for  the  orthotropic  laminates  were  compared  with  the  quasi-isotropic 
laminate  which  served  as  a  reference  case.  Similarly,  the  w/d  and  e/d 
effects  were  evaluated  relative  to  an  infinite-laminate  reference  case. 


SYMBOLS 

d  diameter  of  hole,  m 

^xx'^yy  Young’s  modulus  in  x-  and  y-directions,  respectively,  MPa 
e  edge  distance,  m 

Gyry  shcar  modulus,  MPa 

K3I3  shear  stress-concentration  factor,  (^xy^max/^b 

(^sb^oo  ^sb  infinite  laminate 

tensile  stress-concentration  factor,  (^90)max/^b 
(^tb^co  ^tb  infinite  laminate 

K^g  gross-section  stress-concentration  factor,  (^09)max/Sg 

P  pin  load,  N 

R  radius  of  hole,  m 

r,0  polar  coordinates,  m  and  deg 

Sj3  nominal  bearing  stress,  P/td,  MPa 

Sg  nominal  gross-section  stress,  MPa 

t  laminate  thickness,  m 

w  width  of  laminate,  m 

x,y  Cartesian  coordinates,  m 


2 


0 

9 

V 

a 


c 

m 

xy 


(^90 )  max 
'^xy 

(^xy) max 


one-half  the  angle  of  contact  between  pin  and  hole,  deg 
angular  location  for  (^00) max  hole  boundary,  deg 
Poisson's  ratio 
stress,  MPa 

polar  stress  components,  MPa 

maximum  tangential  stress  on  hole  boundary,  MPa 
Cartesian  shear  stress,  MPa 

maximum  on  shear-out  plane  (y  =  R)  ,  MPa 


FINITE-ELEMENT  ANALYSIS 

Figure  1  shows  the  general  problem  analyzed  in  this  study.  The  laminate 
had  a  frictionless  hole  which  was  loaded  by  a  steel  pin  of  the  same  diameter. 

The  laminate  reacted  to  the  pin  load  by  a  uniform  stress  on  its  left  end.  The 
hole  diameter  was  unity  and  the  various  combinations  of  w/d  and  e/d  were 
generated  by  varying  w  and  e  in  the  finite-element  models.  The  w/d 
effects  were  analyzed  using  a  single,  large  value  of  e/d  (e/d  =  10).  Simi- 
larly,  the  e/d  effects  were  studied  using  w/d  =  20.  For  all  cases,  the  hole 
was  at  a  distance  of  lOd  from  the  left  (loaded)  end  of  the  laminate. 

Table  I  presents  the  elastic  constants  for  the  six  laminates  analyzed. 

These  laminates  were  selected  to  cover  a  broad  range  of  composite  properties. 

The  0°  laminate  properties  (ref.  1)  were  used  with  laminate  theory  to  calculate 
the  elastic  constants  for  the  other  laminates.  Although  the  [ 0^/145^/90^] s 
laminate  was  called  a  quasi- isotropic  laminate,  its  calculated  elastic  constants 
satisfy  all  the  criteria  for  isotropy.  Consequently,  the  results  for  the  quasi¬ 
isotropic  case  were  compared  with  isotropic  solutions  from  the  literature  as 
a  check  on  the  adequacy  of  the  finite-element  analysis. 

The  finite-element  model  for  one-half  of  a  20d  by  20d  square  laminate  is 
shown  in  figure  2.  To  illustrate  the  model  details,  it  is  shown  in  four  parts. 
This  was  the  largest  model  analyzed  (w/d  =  20  and  e/d  =  10).  Smaller  values 
of  w/d  and  e/d  were  obtained  by  removing  elements  from  the  width  and 
unloaded  end  of  the  model,  respectively.  This  model  employed  two  types  of 
finite  elements.  Linear-strain  (6-node)  elements  comprised  the  2.5d  by  5d 
rectangular  region  around  the  hole.  Uniform-strain  (3-node)  elements  comprised 
the  part  of  the  model  outside  this  region.  The  model  was  analyzed  using  the 
BEND  finite-element  computer  program  described  in  reference  2. 

Note  that  the  pin  is  also  modeled  in  figure  2.  The  pin  was  loaded  at 
its  center  and  was  connected  to  the  laminate  by  short,  stiff  spring  elements 
(0.002d  long)  between  adjacent  nodes  across  the  pin/laminate  interface.  Because 
these  springs  had  no  transverse  stiffness,  they  transferred  only  radial  loads 
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and  thereby  produced  the  desired  frictionless  interface*  Also,  these  springs 
were  used  to  allow  the  pin  and  laminate  to  separate  along  part  of  their  inter¬ 
face,  This  was  accomplished  by  assigning  a  zero  stiffness  to  springs  along  the 
separation  arc.  References  3  and  4  show  that  this  separation  V7as  independent 
of  loading  level  for  the  special  case  of  an  initially  snug  fit  (zero  clearance) 
between  the  pin  and  laminate.  For  this  reason,  a  snug  fit  was  used  in  the 
present  study.  For  each  combination  of  material  and  configuration,  the  separa¬ 
tion  angle  was  estimated  and  a  computer  run  was  made.  When  a  spring  was  com¬ 
puted  to  have  a  tensile  force,  its  stiffness  was  set  equal  to  zero  and  that 
case  was  rerun.  Such  reruns  were  needed  for  only  a  few  cases. 


RESULTS  AND  DISCUSSION 
Evaluation  of  Finite-Element  Procedure 

The  finite-element  procedure  was  evaluated  by  comparing  results  for  the 
quasi- isotropic  laminate  with  results  in  the  literature  for  an  isotropic, 
infinite  sheet  containing  a  rigid  pin  with  a  snug  fit.  The  appendix  describes 
the  present  finite-element  analysis  used  to  approximate  this  infinite-sheet  case. 
The  results  for  this  special  case  are  presented  in  figure  3  as  stress  distri¬ 
butions  around  the  hole.  The  stresses  are  normalized  by  the  nominal  bearing 
stress  S]3,  which  is  defined  as  P/td,  where  P  is  the  pin  load,  t  is  the 
laminate  thickness,  and  d  is  the  hole  diameter. 

Figure  3  also  includes  a  curve  for  =  -4/7T  cos  0  because  several 

authors  have  assumed  this  cosine  distribution  for  bearing  stress,  either 
directly  or  indirectly,  by  assuming  a  rigid-body  displacement  for  the  pin.  The 
present  results  (solid  curve)  and  the  results  of  reference  5  (dashed  curve)  are 
in  general  agreement  with  the  cosine  distribution,  except  for  the  angle  of 
contact  between  the  pin  and  the  laminate.  The  90^  contact  angle  9^^  for  the 
cosine  distribution  is  slightly  larger  than  the  83^  contact  angle  found  in  the 
present  study  and  in  references  3  and  5. 

The  000  curves  reach  a  peak  (cr00)iTiax  the  end  of  the  contact  arc. 

This  peak  was  represented  by  a  tensile  stress-concentration  factor  In 

references  6  and  7,  it  was  assumed  that  0^  =  90°,  and  K^b  was  calculated 
as  0.81.  In  references  3  and  5,  however,  K^b  was  calculated  as  0.92  and  0.91, 
respectively,  when  0^  =  83°  was  used.  As  shov/n  in  figure  3,  these  results 
bracket  the  present  K^b  0.87,  In  general,  the  present  results  agree 
reasonably  well  with  the  available  solutions  for  pin-loaded  holes  in  isotropic 
sheets . 

As  an  additional  check  on  the  present  finite-element  model,  stresses  were 
computed  for  an  open  hole  (zero  bolt  modulus)  in  a  sheet  under  remote  uniaxial 
tension.  For  isotropic  material  properties,  a  gross-section  stress-concentration 
factor  K^g  of  3.01  was  calculated,  compared  with  the  familiar  value  of  3.00  for 
an  infinite  sheet.  For  the  extreme  case  of  0°  laminate  properties,  K^g  =  6.45, 
which  compares  favorably  with  the  infinite-sheet  value  of  6.43  from  reference  8. 
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width  Effects 


As  previously  mentioned,  six  laminates  were  analyzed.  The  w/d  values 
ranged  from  2  to  20  and  e/d  was  10.  For  each  case,  the  laminate  was  loaded 
by  a  steel  pin  having  a  snug  fit.  Results  are  first  presented  as  stress  dis¬ 
tributions  near  the  hole  and  are  briefly  discussed  in  terms  of  the  basic  fail¬ 
ure  modes:  net-tension,  shear-out,  and  bearing.  Next,  stress-concentration 
factors  are  plotted  for  a  range  of  w/d  values.  Finally,  these  stress- 
concentration  factors  are  compared  with  their  corresponding  inf inite- laminate 
values.  The  quasi-isotropic  laminate  serves  as  a  reference  case  to  evaluate 
the  effects  of  anisotropy. 

Stress  distributions.-  Stresses  for  the  quasi-isotropic  laminate  are  shown 
in  figure  4.  For  simplicity,  only  results  for  w/d  =2,  3,  and  oo  are  shown; 
other  results  for  the  full  w/d  range  are  summarized  in  table  II.  In  fig¬ 
ure  4(a),  the  ^199  curve  for  w/d  =  2  has  a  maximum  value  (^00) max 
more  than  twice  that  for  w/d  =  In  contrast,  the  curves  show  only 

minor  differences  for  the  various  w/d  values.  In  fact,  the  maximum  is 

nearly  unaltered,  and  the  contact  angle  9^  varies  only  from  about  83^  to  90^, 

Figure  4(b)  shows  the  shear  stress  computed  for  the  shear-out  failure 

plane,  y  =  R.  The  narrow  laminate  width  of  w/d  =  2  elevated  this  T^y  by 
about  25  percent  above  the  infinite-laminate  case. 

The  ^QQ,  ^rr'  ^xy  figure  4  are  associated,  respectively,  with 

the  following  failure  modes:  net-tension,  shear-out,  and  bearing.  Of  these, 
only  cr99  shows  much  dependence  on  w/d.  As  a  result,  only  the  net-tension 
failure  strength  should  be  strongly  influenced  by  laminate  width. 

The  0^  laminate  stresses  are  presented  in  figure  5.  The  a99  curves  for 
this  highly  orthotropic  case  have  larger  maximum  values  and  larger  gradients 
than  for  the  quasi-isotropic  case.  However,  (^90) max  ^  only 

about  twice  as  large  as  the  corresponding  value  for  w/d  =  <»;  this  w/d  effect 

is  comparable  to  that  found  for  the  quasi- isotropic  laminate.  The  Oj.j.  curves 

differ  markedly  from  the  cosine  distribution  and  have  a  sharp  peak  at  0  =  0*^, 
which  is  the  stiff est  direction  for  this  0^  laminate.  The  curves  have 

almost  no  dependence  on  w/d  and  show  a  contact  angle  0^  of  about  82*^.  The 
x^y  curves  are  also  only  slightly  different  for  the  range  of  w/d  values. 

Results  for  the  90^  laminate  are  presented  in  figure  6.  Although  the  099 
curves  have  the  expected  peaks  near  0  =  90°,  they  have  higher  peaks  at  0  =  0°. 
Note  that  in  this  90°  laminate  the  fibers  are  tangent  to  the  hole  at  the  point 
of  maximum  CJ99 .  Again,  the  curves  are  quite  similar  to  one  another;  all 

deviate  from  the  cosine  distribution.  Also,  they  deviate  significantly  from  the 
two  previous  laminates  by  having  their  maxima  near  0  =  70°.  The  contact  angles 
in  figure  6  range  from  about  84°  to  86°.  The  t^y  curves  in  figure  6  have 
higher  peaks  than  those  found  for  the  two  previous  laminates. 

The  [0°/90°]s  laminate  results  are  plotted  in  figure  7.  As  expected,  these 
results  are  approximately  equal  to  averages  of  those  shown  for  0°  and  90°  lami¬ 
nates.  The  (^99) max  value  for  w/d  =2  is  located  at  0  =  90°  and  is  about 
twice  as  large  as  its  infinite-laminate  value.  The  curves  and  the  x^y 

curves  are  like  those  for  the  quasi- isotropic  case  in  figure  4. 
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The  stress  distributions  for  the  [ ±45]  g  laminate,  plotted  in  figure  8,  are 
unlike  those  in  previous  figures.  The  099  curve  for  w/d  =  «>  has  a  pro¬ 
nounced  peak  at  0  =  45°,  where  the  -45°  ply  fibers  are  tangent  to  the  hole. 

For  smaller  w/d  values,  the  CTqq  curves  have  a  second  peak  near  0  =  135° 
where  the  45°  ply  fibers  are  tangent  to  the  hole.  This  second  peak  appears  to 
be  strongly  amplified  for  narrow  laminates.  Furthermore,  narrow  laminates  show 
larger  contact  angles  and  larger  stresses  near  the  end  of  the  contact  arc  than 
wide  laminates. 

The  [o°/±45°]g  results  shown  in  figure  9  have  the  expected  trends.  For 
small  w/d  values,  (°00)max  nearly  three  times  the  infinite-laminate  value. 
This  peak  value  occurs  slightly  beyond  the  region  of  contact.  The  Oj-j-  and 
Txy  curves  are  similar  to  the  quasi- isotropic  results  and  are  not  strongly 
influenced  by  the  laminate  width. 

Stress-concentration  factors.-  To  compare  the  six  different  laminates,  the 
stress-concentration  factor  Kj-b  was  calculated  using  each  value  of  (°09)raax 
in  the  previous  figures.  These  K^b  values  are  listed  in  table  II  along  with 
the  location  0jjj  for  each  value  of  (a0Q)jj|gjj  and  the  contact  angle  0,^  that 
occurred  for  each  case.  The  K^b  values  are  also  plotted  in  figure  10  for  the 
range  2  i  w/d  i  10.  This  figure  shows  similar  width  effects  for  all  six  lami¬ 
nates;  the  larger  w/d  values  produced  smaller  K^-b  values.  As  previously 
mentioned,  the  quasi-isotropic  case  provides  a  reference  for  evaluating  the 
anisotropy  effects.  The  curve  for  this  reference  case  lies  near  the  lower 
extreme  for  the  range  of  computed  results.  Only  the  [ 0°/±45°]  g  curve  is  lower. 
Note  that  the  results  from  reference  9  for  a  finite— width  isotropic  sheet  agree 
closely  with  the  present  quasi-isotropic  results.  As  expected,  the  0°  laminate 
curve  is  the  upper  bound  for  the  range  of  results,  and  the  [ 0°/90°] g  curve  lies 
near  the  middle  of  this  range.  Comparisons  involving  the  90°  laminate  and  the 
[  ±45°]  g  laminate  are  difficult  because  their  (a09)jj,3j^  values  occurred  at 
9n,  =  0°  and  45°,  respectively,  rather  than  at  the  usual  0;^  «  90°  for  the 
laminates  containing  some  0°  plies. 

The  curves  in  figure  10  appear  to  have  a  negative  slope  even  at  w/d  =  10, 
suggesting  that  the  width  effects  may  exist  beyond  w/d  =10.  To  investigate 
this,  the  Kbb  values  were  each  normalized  by  their  corresponding  infinite- 
laminate  value  (K(.b)oo  and  then  replotted  in  figure  11.  Of  the  six  curves  in 

figure  11,  only  those  for  the  90°  laminate  and  the  [ ±45°] 3  laminate  approach 
unity,  indicating  that  only  for  these  laminates  is  the  effect  negligible  for 
w/d  =  10.  The  other  four  curves  converge  to  about  1.17.  This  means  that  K^b 
is  amplified  by  17  percent  even  for  a  rather  wide  plate  with  w/d  =  10.  For 
comparison,  if  the  present  laminates  had  unloaded  holes  and  remote  uniaxial 
loads,  their  stress  concentrations  would  be  amplified  by  less  than  about  2  per¬ 
cent  at  w/d  =  10  (ref.  10).  This  comparison  shows  that  a  finite  width  influ¬ 
ences  the  tensile  stress  concentration  at  a  loaded  hole  much  more  than  it  does 
for  an  unloaded  hole.  In  figure  11,  the  curves  for  all  the  laminates  that  con¬ 
tain  some  0°  plies  lie  close  to  the  quasi- isotropic  reference  curve.  In  fact, 
for  w/d  ^  3,  use  of  the  quasi-isotropic  curve  to  approximate  the  others  with 
0°  plies  would  result  in  errors  of  less  than  10  percent.  The  other  two  lami¬ 
nates,  90°  and  [±45°]s,  have  curves  in  figure  11  that  lie  well  below  that  for 
the  quasi- isotropic  case.  However,  these  two  laminates  have  their  peak  values 
at  0j„  =  0°  and  45°,  respectively. 
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Edge-Distance  Effects 


As  previously  described,  the  six  laminates  were  analyzed  for  edge-distance 
effects  by  using  a  value  of  w/d  =  20.  The  e/d  values  ranged  from  1  to  20. 

The  laminates  were  assumed  to  be  loaded  by  a  steel  pin  fitting  snugly  in  the 
hole.  First,  the  results  are  plotted  as  stress  distributions  for  each  laminate. 
The  stress-concentration  factors  for  the  six  laminates  are  then  compared  with 
one  another  and  with  their  corresponding  infinite-laminate  values. 

Stress  distributions.-  The  quasi- isotropic  case  was  again  used  as  a  refer¬ 
ence  to  investigate  the  effects  of  laminate  anisotropy  on  stresses  at  a  loaded 
hole.  Stress  distributions  for  the  quasi- isotropic  case  are  presented  in  fig¬ 
ure  12.  The  curves  for  two  finite  edge  distances  (e/d  =  1  and  2)  are  shown 
along  with  those  for  an  infinite  laminate  (e/d  =  c»)  ,  For  the  extreme  e/d  =  1 
case  in  figure  12(a),  the  peak  aQQ  is  about  twice  as  large  as  that  for 
e/d  =  00.  As  expected,  this  peak  occurs  near  the  end  of  the  contact  arc.  But 
for  small  e/d  values,  this  peak  shifts  slightly  to  the  left,  in  contrast  to 
figure  4(a),  which  shows  a  slight  shift  to  the  right  for  small  w/d  values. 
Otherwise,  however,  the  099  trends  in  figure  12(a)  are  quite  similar  to  those 
in  figure  4(a).  The  Oj-j-  and  curves  in  figures  12(a)  and  12(b)  show 

little  influence  of  e/d,  as  long  as  e/d  ^  2. 

Figure  13  presents  stresses  computed  for  the  0^  laminate.  This  extremely 
orthotropic  laminate  is  strongly  influenced  by  finite  edge  distances.  Even  for 
e/d  =  2,  figure  13(a)  shows  that  (^00)niax  about  twice  as  large  as  for  the 
reference  case  (e/d  =  00).  For  e/d  =  1,  it  is  about  three  times  as  large. 

(See  table  III.)  This  stress  amplification  due  to  finite  e/d  values  is  larger 

than  that  shown  in  the  previous  section  for  finite  w/d  values.  This  is  prob¬ 

ably  attributable  to  the  fact  that  the  0^  laminate  is  very  stiff  in  the  edge- 
distance  direction  but  very  flexible  in  the  width  direction.  In  contrast  to 
the  cr99  curves,  the  curves  show  very  little  e/d  influence;  the  contact 

angle  0^  is  about  82^  for  the  full  range  of  e/d  values. 

Figure  13(b)  shows  that  e/d  has  a  strong  influence  on  the  T^y  values 
computed  for  the  shear-out  plane  (y  =  R) .  The  curve  for  e/d  =  1  has 

a  peak  that  is  about  three  times  as  large  as  the  inf inite- laminate  value.  The 
(Txy)max  values  are  similar  for  both  the  0^  and  quasi-isotropic  laminates.  Of 
course,  these  laminates  have  different  shear  strengths  and  therefore  should 
respond  differently  to  the  same  T™  stress.  For  example,  the  0^  laminate 
has  a  shear  strength  that  is  only  about  one-fourth  that  for  the  quasi- isotropic 
case  and  is  therefore  more  vulnerable  to  the  shear-out  failure  mode. 

The  90^  laminate  stresses  are  plotted  in  figure  14.  All  three  stresses 
show  relatively  little  influence  of  e/d.  These  stresses  were  increased  by  only 
about  50  percent  even  for  the  extreme  e/d  =  1  case.  Perhaps  this  laminate 
shows  little  sensitivity  to  the  free  edge  because  it  has  a  very  low  stiffness 
in  the  edge-distance  (x-axis)  direction. 

The  [O^/SO^lg  laminate  is  relatively  stiff  in  the  x-direction  and,  as 
expected,  has  a  much  stronger  dependence  on  e/d  than  observed  for  the  90*^ 
laminate.  For  e/d  =  1  (fig.  15),  the  (a99)max  is  about  twice  as  large 
as  when  e/d  =  00.  Also,  {cr99)niax  occurs  at  0  =  0°  for  this  e/d  =  00 
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reference  case  but  occurs  at  6  =  90°  for  finite  edge  distances.  Similarly, 
the  Gj-r  peak  occurs  at  6=0°  for  e/d  =  co  but  shifts  toward  9  =  45°  for 
finite  e/d  values. 

Figures  16  and  17  present  stress  distributions  for  the  [+45°  ]g  and  the 
[0°/±45°]s  laminates,  respectively.  As  expected,  the  f±45°]s  curves  in  fig¬ 
ure  16(a)  have  their  peak  value  at  0  =  45°,  where  the  -45°  ply  fibers  are 
tangent  to  the  hole.  This  highly  stressed  -45°  ply  is  probably  responsible  for 
the  relatively  large  T™  peak  shown  in  figure  16(b).  The  Oqq  curves  for 
the  [0°/±45°]3  laminate  have  their  peak  at  the  usual  0  =  90°  location.  In 
fact,  all  the  [0°/±45°]g  results  in  figure  17  are  quite  similar  to  the 
[(D°/±45°/90° ]g  reference  case  in  figure  12.  Apparently,  the  90°  plies  have 
little  influence  on  the  stresses  near  a  loaded  hole  when  they  are  used  together 
with  0°  and  ±45°  plies. 

Stress-concentration  factors.-  Stress-concentration  factors  for  the 

six  laminates  are  compared  for  a  range  of  e/d  in  figure  18.  These  K^b  values 
were  calculated  using  (°09)inax  t^°"*  °60  stress  distributions.  These  Kj-b 

values  are  also  presented  in  table  III,  along  with  the  contact  angle  and 

the  angular  location  0]j|  for  each  case.  Figure  18  shows  that,  despite  the  wide 
range  of  K^b  for  e/d  =  1,  all  six  curves  lie  within  a  range  of  0.90  to  2.16 
for  e/d  =  10.  As  in  the  previous  analysis  for  w/d  effects,  the  quasi¬ 
isotropic  reference  curve  lies  near  the  lower  end  of  this  range.  Apparently, 
the  quasi- isotropic  laminate  provides  the  most  uniform  load  transfer  near  a 
loaded  hole.  Anisotropy  produces  more  concentrated  load  paths  and  correspond¬ 
ingly  higher  stress  concentrations. 

Figure  19  presents  a  comparison  of  these  K^-b  results  in  a  different  way. 
Each  Ktb  value  was  normalized  by  the  corresponding  (Kt.b)co  value  for  that 
laminate.  Each  curve  in  figure  19  should  approach  unity  for  sufficiently  large 
e/d  values,  and  they  all  appear  to  have  this  trend.  For  e/d  =  10,  each  curve 
lies  within  the  narrow  range  of  1.0  to  1.1.  However,  for  e/d  =  3,  a  practical 
value  for  bolted  joints,  the  Kj-j.,  for  a  quasi-isotropic  laminate  is  about 
25  percent  larger  than  for  the  infinite-laminate  case;  Ktb  values  for  all  the 
other  laminates  with  0°  plies  are  even  larger.  The  largest  increase  is  about 
50  percent  for  the  0°  laminate. 

Because  shear-out  is  a  viable  failure  mode  for  finite  e/d  cases,  stress- 
concentration  factors  were  also  calculated  for  the  stress  along  the  shear- 

out  planes.  Figure  20  shows  the  stress-concentration  factor  Kgb  =  ('rxy)max/Sb’ 
(Also  see  table  IV.)  For  e/d  =  1,  the  Kgb  results  are  rather  tightly  grouped 
in  the  range  0.7  to  0.9.  But  for  larger  e/d  values,  Kgb  ranges  from  about 
0.2  to  0.6.  As  expected,  the  90°  laminate  provides  the  upper  bound  for  this 
range  and  the  0°  laminate  provides  the  lower  bound.  The  quasi- isotropic  curve 
lies  near  midrange. 

Each  Kgb  was  normalized  by  the  corresponding  infinite- laminate  value 
(^sb^oo*  These  Kgb/(Kgb)oo  results  are  shown  in  figure  21.  As  expected,  for 
small  e/d  values,  the  results  have  a  wide  range.  But  for  e/d  S  3,  all  six 
curves  are  closely  grouped  and  approach  unity.  For  e/d  >  3,  the  e/d  effect 
is  less  than  about  25  percent,  even  for  the  extremely  orthotropic  0°  laminate. 
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Stress-concentration  factors  for  the  a^-j-  stress  were  also  calculated, 
but  the  e/d  effects  were  found  to  be  quite  small.  For  e/d  ^  3,  the  peak 
stresses  were  less  than  10  percent  larger  than  corresponding  infinite- 
laminate  values.  Consequently,  e/d  should  have  a  negligible  effect  on  bear¬ 
ing  strength  for  e/d  ^  3. 


CONCLUDING  REMARKS 

Stresses  were  calculated  for  finite-size  orthotropic  laminates  that  are 
loaded  by  a  frictionless  steel  pin  in  a  circular  hole  of  the  same  diameter. 

These  calculations  were  based  on  finite-element  analyses  for  six  laminates: 
quasi- isotropic  [  0^/145^/90^]  sf  90^,  [  0^/90^]  sr  [l45^]s/  and  [0^^/1450]  g, 

Stress-concentration  factors  were  based  on  nominal  bearing  stress  and  were 
determined  for  wide  ranges  of  width- to-diame ter  ratio  w/d  and  edge-distance- 
to-diameter  ratio  e/d.  For  reference  purposes,  the  infinite-laminate  case 
was  also  analyzed  for  each  laminate. 

When  load  is  introduced  into  a  laminate  by  a  pin,  the  hole  surface  and 
pin  remain  in  contact  only  over  a  portion  of  the  hole  boundary  defined  by  a 
contact  arc.  In  general,  for  cases  with  either  an  initial  clearance  fit  or  an 
interference  fit,  this  contact  arc  varies  with  load.  However,  recent  liter¬ 
ature  shows  that  the  present  case  with  an  initially  snug  fit  has  a  contact  arc 
that  is  independent  of  load.  Although  each  combination  of  laminate  material 
and  configuration  produced  a  different  contact  arc,  each  such  arc  was  in  the 
narrow  range  of  160°  to  180°.  Therefore,  neither  the  wide  ranges  of  material 
anisotropy  nor  laminate  configuration  had  much  influence  on  the  contact  arc. 

As  expected,  anisotropy  had  a  strong  influence  on  the  contact  stress  distri¬ 
butions.  These  distributions  deviated  significantly  from  the  well-known  cosine 
distribution,  especially  for  the  0^  and  90^  laminates. 

The  anisotropy  also  had  a  strong  influence  on  the  magnitude  and  location 
of  the  peak  tangential  stress  on  the  hole  boundary.  This  peak  was  represented 
by  a  stress-concentration  factor  K^j^.  For  the  infinite-laminate  case,  the 
stress  concentration  factor  ranged  from  0.82  to  2.16  compared  to  0.87  for  the 
quasi-isotropic  laminate.  The  stress  peak  developed  near  the  end  of  the  contact 
arc  except  for  the  90°  and  [ ±45°] g  laminates,  which  had  their  peaks  where  the 
ply  directions  (fibers)  were  tangent  to  the  hole. 

The  narrow  laminates  produced  large  values  for  the  loaded  holes.  Even 

for  a  rather  wide  (w/d  =  10)  configuration,  the  quasi-isotropic  laminate  had  a 

that  was  nearly  20  percent  larger  than  its  infinite-laminate  reference  value. 
In  contrast,  for  an  unloaded  hole,  the  corresponding  stress-concentration  factor 
is  only  about  2  percent  larger  than  its  infinite-laminate  value.  In  general, 
the  orthotropic  laminates  containing  0^  plies  produced  trends  quite  like  that  for 
the  quasi- isotropic  laminate,  especially  for  w/d  ^  3.  Also  for  w/d  ^  3,  the 
peak  bearing  stresses  for  all  six  laminates  were  elevated  by  less  than  about 
10  percent. 

As  expected,  the  finite-edge  distance  also  increased  the  stresses  near  the 
loaded  hole.  However,  this  influence  varied  widely  for  the  different  laminates 
and  different  stresses.  For  the  extreme  case  of  e/d  =  1,  this  influence 
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increased  the  0^  laminate  K^b  by  more  than  200  percent  but  increased  the 
90^  laminate  K^b  by  only  50  percent.  The  edge  distances  only  weakly 
influenced  the  peak  values  for  shear-out  and  bearing  stress.  For  e/d  ^  3, 
the  peak  shear  stresses  were  increased  by  less  than  about  25  percent.  Simi¬ 
larly,  for  e/d  ^  3,  the  peak  bearing  stresses  were  increased  by  less  than 
10  percent. 

In  summary,  both  finite  width  and  finite  edge  distance  had  a  strong  influ¬ 
ence  on  the  maximum  tensile  stress,  but  had  much  less  influence  on  the  other 
peak  stresses  near  the  hole. 


Langley  Research  Center 

National  Aeronautics  and  Space  Administration 
Hampton,  VA  23665 
April  1,  1981 
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APPENDIX 


APPROXIMATION  FOR  INFINITE- LAMINATE  CASE 

To  analyze  the  infinite-laminate  case,  the  finite-element  procedure  was 
systematically  applied  for  increasing  values  of  w/d.  The  computed  stress- 
concentration  factor  Ktb  for  a  loaded  hole  was  expected  to  asymptotically 
approach  an  infinite-laminate  value  (Ktb)oo-  ancillary  "convergence” 

study  was  conducted  with  a  square,  quasi-isotropic  laminate  and,  for  com¬ 
parison  with  the  literature,  the  pin  was  assumed  to  be  rigid. 

Figure  22  and  table  V  present  K^b  values  for  the  three  different  load¬ 
ing  cases  analyzed  in  the  convergence  study.  Case  1 ,  represented  by  the  upper 
curve,  is  the  basic  problem  dealt  with  in  this  study.  Case  2  is  included  for 
comparison  purposes.  Note  that  cases  1  and  2  have  the  remote  stress  reactions 
on  opposite  ends  of  the  model.  Otherwise,  these  two  cases  are  identical  and, 
for  an  infinite  laminate,  they  should  have  the  same  (Ktb)oo  value  because  the 
remote  stress  would  be  zero  for  both  cases.  On  the  other  hand,  the  difference 
between  the  Ktb  curves  for  cases  1  and  2  is  a  measure  of  the  finite-size 
effect  for  a  given  w/d  value.  The  curves  for  cases  1  and  2  converge  slowly 
in  figure  22  and  differ  by  almost  20  percent  at  w/d  =  20. 

Because  neither  the  case  1  curve  nor  the  case  2  curve  approached  a  limit 
for  the  present  range  of  w/d,  the  third  case  was  considered.  This  case  could 
be  generated  by  superimposing  cases  1  and  2.  Therefore,  if  cases  1  and  2  have 
the  same  limit,  then  the  limit  must  also  apply  to  case  3.  Conversely,  a  limit 
found  for  case  3  would  also  be  the  desired  limit  for  case  1.  Toward  this  end, 
case  3  loading  was  used  in  the  finite-element  analysis  for  a  range  of  w/d 
values.  The  resulting  K^b  curve  in  figure  22  has  nearly  converged  for 
w/d  =  20. 

To  check  the  assumption  that  all  three  curves  in  figure  22  converge  to 
(Ktb)oo^  each  curve  was  extrapolated  by  using  a  numerical  analysis.  Each  curve 
was  approximated  by  the  cubic  polynomial 


Ktb  -  ^0 


The  coefficients  aQ  to  a3  were  obtained  by  a  least-squares  fit  to  the 
values  in  table  V.  The  polynomials  were  each  found  to  have  a  correlation  coef¬ 
ficient  that  deviated  from  unity  by  less  than  4  x  10“^,  indicating  nearly 
perfect  agreement  with  the  tabulated  values.  For  very  large  w/d  values, 

each  polynomial  approaches  its  constant  ag  term.  Therefore,  each  ag  was 
an  estimate  for  (Ktb)cx>*  three  aQ  values  were  0.863,  0.864,  and  0.863, 

respectively,  for  cases  1,  2,  and  3.  Therefore,  the  three  curves  converge  to 
the  same  limit  of  about  0.863. 
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For  convenience,  the  case  3  results  for  w/d  =  20  were  assumed  to 
approximate  the  infinite-laminate  results  for  case  1.  The  resulting  value 
(Ktb)  oo  =  0.869  was  considered  to  be  sufficiently  close  to  the  0.863  extrapolated 
limit.  The  infinite-laminate  stress  distributions  and  ^  values  for  the 

other  laminates  were  also  approximated  using  case  3  loading  with  w/d  =  20. 

The  (Ktb)oo  values  are  presented  in  table  II  and,  for  comparison  purposes,  are 
repeated  in  table  III.  The  shear  stress-concentration  factors  (Ksb)  ^or  an 
infinite  laminate  are  presented  in  table  IV. 
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TABLE  I.-  LAMINATE  ELASTIC  CONSTANTS 


Laminate 

Elastic  constants 

Exx'  MPa 

Eyyf  MPa 

Gjjy,  MPa 

Quasi- isotropic 
[0O/±45°/90°]s 

57  890 

57  890 

22  090 

0.310 

0° 

146  900 

10  890 

6  412 

0.380 

90° 

10  890 

146  900 

6  412 

0.028 

[0°/90°]s 

79  500 

79  500 

6  412 

0.052 

[±45°]s 

22  250 

22  250 

37  770 

0.735 

[0°/±45°]s 

64  200 

28  210 

27  320 

0.699 

TABLE  II.-  TENSION  STRESS-CONCENTRATION  FACTORS  FOR  w/d  RANGE  WITH  e/d 
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TABLE  IV.-  SHEAR  STRESS-CONCENTRATION  FACTORS  FOR  SHEAR-OUT  PLANE  WITH  e/d  RANGE  AND  w/d 
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TABLE  V.-  STRESS-CONCENTRATION  FACTORS  FOR  A  SQLJARE  QUASI- ISOTROPIC 


LAMINATE  LOADED  BY  A  RIGID  PIN;  9^  =  83^ 


Remote  loading  type^ 


Tension 
(case  1) 


Canpression 
(case  2) 


Tension/ compression 
(case  3) 


K 


tb 


<3eg 


Ktb 


6m,  deg 


K 


tb 


deg 


1.832 

^(1.728) 


90 


0.576 

^(0.542) 


82 


1 .157 
^(1.144) 


90 


1.304 

^(1.283) 


90 


0.693 

^(0.645) 


82 


0.964 


86 


1.040 


86 


0.752 

^(0.736) 


82 


0.888 


86 


0.975 


86 


0.777 

^(0.774) 


82 


0,874 


86 


0.944 


86 


0.794 


86 


0.869 


86 


^See  figure  22  for  description  of  loading  types, 
^Stress  concentration  at  0  =  86*^. 


0,  deg 


Figure  3.-  Stresses  in  quasi-isotropic  laminate  loaded  by  a  rigid  pin 


(a)  Stresses  along  hole  boundary  (r  =  R) . 


(b)  Shear  stress  along  shear-out  plane  (y  =  R) . 

Figure  4,-  Stresses  in  quasi-isotropic  laminate  as  a  function  of 

w/d;  e/d  =  10. 


(a)  Stresses  along  hole  boundary  (r  =  R) - 


(b)  Shear  stress  along  shear-out  plane  (y  =  R) . 

Figure  7.-  Stresses  in  fo^/90^]g  laminate  as  a  function  of  w/d;  e/d 


(a)  Stresses  along  hole  boundary  (r  =  R) . 


Figure  9.-  Stresses  in  [ 0^/45^] s  laminate  as  a  function  of  w/d 


Figure  14.-  Stresses  in  90^  laminate  as  a  function  of  e/d;  w/d 


(a)  Stresses  along  hole  boundary. 


Figure  17.-  Stresses  in  [0^/±A5^] ^  laminate  as  a  function  of  e/d;  w/d 


Figure  18,-  uqq  stress-concentration  factors  as  a  function  of 

e/d;  w/d  =  20. 


0  2  4  6  8  10 

e/d 


Figure  19.  Normalized  CTqq  stress-concentration  factors  as  a  function  of 

e/d;  w/d  =  20. 


0  2  4  6  8  1( 

e  /d 


Figure  21 Normalized  x^y  stress-concentration  factors  for  shear-out 

plane  (y  =  R) ;  w/d  =  20. 
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